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The general concern of this paper is the effect of rough boundaries on fluids. We consider 
a stationary flow, governed by incompressible Navier-Stokes equations, in an infinite domain 
bounded by two horizontal rough plates. The roughness is modeled by a spatially homoge- 
neous random field, with characteristic size e. A mathematical analysis of the flow for small 
£ is performed. The Navier's wall law is rigorously deduced from this analysis. This extends 
substantially former results obtained in the case of periodic roughness, notably in |151 116j . 
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1 Introduction 

The understanding of roughness-induced effects is a major concern in fluid dynamics. Indeed, 
many examples of physical relevance involve rough boundaries. By "rough", we mean that 
the spatial variations are small compared to the typical length of the problem. For instance, 
in geophysics, the bottom of the oceans and the shores are rough with respect to the large 
scale flow. Also, in an industrial framework, containers have often imperfections that qualify 
them as rough. 

The main problem is to know in which way such irregular boundaries affect the flow. 
This is especially important with regards to numerical computations: indeed, roughness is in 
general too small to be captured by the discretization grid of the simulations. 

To overcome this difficulty, one often relies on wall laws. A wall law is a boundary 
condition that is imposed on an artificial boundary inside the domain. The idea is to filter 
out the precise description of the flow near the real rough boundary. The wall law should 
only reflect the large scale effect of the roughness, in the spirit of a homogeneization process. 

In many cases, the determination of wall laws relies on formal calculations, grounded by 
empirical arguments (see for instance jHl HOI)- The present paper is a mathematical justi- 
fication of some wall laws, in the case of an incompressible viscous fluid. We consider the 
two-dimensional stationary Navier-Stokes equations: 



where the lower and upper boundaries jf and 7^ are to be precised. As usual, the flelds 
u = {ui,U2){x) £ M?, p = p{x) G M are the velocity and the pressure, and 1/ > is the 
kinematic viscosity. Equations (|l.lj) are completed with the classical no-slip conditions 



M • Vu + - z^An = 0, X G r2^, 
div u = 0, X ^Vf , 




in a domain of channel type: 



= {(xi,X2) G M^ 7f(xi) <X2< ll{xi)} 



(1.2) 



Moreover, we prescribe the fluid flux through the channel, that is condition 




(1.3) 



where </> > is a given constant flux, and 



'^{^i) = Hixi) <X2< ll{xi)] 
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is a vertical section of the channel at xi. Remark that by incompressibility and boundary 
condition ()1.2|) . the left hand-side of does not depend on xi. The functions 7f and 7^ 
model rough plates, with small characteristic size e. Broadly, they read: 

7f = -^li {xi/e), fu = l + e-iu (xi /e) , 

for Lipschitz functions 

11 = iiivi) e]0, l[, 7„ = 7n(yi) e]0, l[. 

More precise assumptions on 7^ and 7^ will be made further on. The set = Mx]0, 1[ will be 
called the interior domain. 

We wish to study solutions of (|1.2() . and to determine appropriate wall 

laws for this system. In other words, we look for operators B^{x,Dx) such that solutions 
of the interior system 

' V -Vv + Vq - v/\v = 0, X G r2, 
div ti = 0, a; G ri, 

t>i = 0, 

<t(xi) 

=0, 

approximate well in J7, for small e. 

The mathematical treatment of wall laws has been the matter of many articles. The note 

is devoted to the analysis of Laplace equation in an annular domain with perforations. 
Numerical and formal computations for fluid flows can be found in j2j An analysis of 
Couette flows in rough domains has been performed in [Sj . Let us also mention the important 
contributions of Jager and Mikelic on wall laws for channel flows ( jl5l fI6j ) . We refer to jl7j on 
a related problem with porous boundaries. Finally, see Jl], [H] for study of roughness-induced 
effects on some geophysical systems. 

All these articles are devoted to periodic roughness, meaning that the boundary functions 
7i and 7„ are periodic. This is of course a mathematical simplification, which is highly 
unrealistic from the point of view of physics. Our goal here is to drop this restriction, and 
treat non-periodic roughness. Precisely, we consider roughness that is distributed following 
a spatially homogeneous random field. A complete description of the rough domain will be 
given in the next section. 

Following |15j in the periodic case, special attention is paid to the simple Dirichlet wall 
law: 

B'{x,D,){v)\Q^ = v\9n = ^, (1.5) 
and to the Navier's friction law: 

B'{x,D^){v)\an = {Ce{x)^-Vr\ \an = 0. (1.6) 



introduced by Navier |2T] and extensively used in simulations of geophysical flows. Losely, 
we show two main results: 



1. The Dirichlet wall law yields a 0(e) approximation of the real solution u^, that is — v' 
is 0{e) in an appropriate quadratic norm, to be described in the next section. 
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2. For appropriate Cs, the Navier's law yields a o(e) approximation of the real solution u'^. 

These results extend those of ^H]- They are deduced for a precise description of for small e, 
especially of the boundary layer flow near dfl^. Precise statements, including the expression 
of Cs, will be given in the next section. 

To end this introduction, let us point out some difficulties related to the proof of these 
results. First, we consider a domain 0^ that is not bounded in the tangential direction 
{xi E M). To our knowledge, previous studies dealt with bounded channel domains, wether 
with lateral boundaries (plus in- and out-flux lateral boundary conditions, see or with 
periodic boundary conditions, see j^l- Note that such periodicity condition is not compatible 
with our non-periodic roughness. Due to the unbounded channel domain, we work with 
only locally integrable functions, which leads to completely different treatment of the energy 
estimates. Secondly, as the roughness is non-periodic, the boundary layer system is more 
complex. Due to the lack of compactness both in the tangential and transverse variables, 
we are not able to solve it in a deterministic setting. We use a variational formulation that 
involves the random variable. In addition to this problem, the behaviour of the boundary 
layer profile far from the boundary is not obvious. This can be understood using formally the 
tangential Fourier transform. Indeed, in the periodic setting, the Fourier modes are discrete, 
and allow a clear separation between the non-oscillating part (the constant mode) and the 
oscillating ones (the non-constant modes). But in the non-periodic case, there is no such 
separation, and Fourier modes close to zero create trouble. To control these low frequencies, 
we must again inject some probabilistic information (namely, the ergodic theorem). This 
difficulty appears also further on in the study, to establish energy estimates. 

The rest of the paper is structured as follows. The next section contains a precise modeling 
of the domain, and the statements of the mathematical results. The third section is devoted 
to the Dirichlet wall law. The fourth section focuses on the boundary layer analysis. The 
final section is the justification of Navier's wall law. 

2 Statement of the results 
2.1 Modeling of the rough domain 

Let e > 0, and {M,M,fi) a probability space. For all m £ M, we define a rough domain 
ri^(m) by 

n%m) = nuRf{m)URi{m), 

where = Mx]0, 1[ is the interior domain, and Rf^{m) is the lower, upper rough part. To 
obtain a realistic model for roughness, we use spatially homogeneous random fields: following 
|22j . |18j or IS], we recall that a homogeneous random field is a measurable map 

7 : M X M" M™ 

satisfying: for all h, zi, . . . , Zj. and all Borel subsets Bi, . . . ,Bj. of M™, 

H {{m G M, 7(m, zi + h) e Bi, . . . , -f{m, Zk + h) e Bt}) 

= /i ({m G M, 7(m, zi) £ Bi,..., 7(771, Zk) £ Bk}) . 
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We remind that for n = m = 1, a homogeneous random field is often called a stationary 
random process. We thus define 

Ri{m) = {x = {xi,X2), xi G M^, > 2:2 > —eji{m,xi/e)] , 
Ru{m) = [x = (3;i,X2), xi G M^, < X2 - 1 < e^u{m,xi/e)] , 

where (7i,7u) = (7«i 7u)("i, yi) £]0, Ip is a homogeneous random field. Moreover, we assume 
that for all m, 7/,u(m, •) is a K-Lipschitz function, with K > independent of m, and that 
^ ^ H,u{^, ■) is measurable with values in the set C6(M; M^) of continuous bounded functions. 

Following a classical construction of Doob (see ^2] or [HI for all necessary details), we 
introduce another probability space, which will be more convenient to our description. Let P 
the set of K-Lipschitz functions uj : M i-^]0, Let V the a-algebra generated by the sets 

{uj e P, cu(yi) e A, Vyi G 5} , 

where i? is a finite subset of Q, and A is a disk with rational center and radius. Note that 
V is simply the borelian o"— algebra of P, seen as a subset of Cb(M; M^). Finally, consider the 
set function vr : 1-^ M given by 

tt{Q) = fi{{me M, (7„,7/)(m,-) G Q}) . 

One can show (c/ ^2]) that vr is a probability measure on {P,V). Moreover, we can define a 
translation group 

Th-P^P, Th{io){yi) = uj{yi + h), Vyi,/iGM, 
that preserves it. In this way, one can also describe the boundaries with the measurable map 

Indeed, the laws of the random variables 

{{7i,lu){-,zi + h),..., {-fi,'yu){-,Zk + h)) 

and 

{{uJi,u;u){-,zi + /i), . . . , {u;i,LLiu)i-,Zk + h)) 

are the same (and independent of h, for any zi, . . . , z^). The advantage of this last framework 
is that one can write 

^u,l{yi,^) = K,l°Ty^{^), with {hu,hi){uj) =uj{0), 

which will be useful in the study of boundary layers. Hence, we will rather use the formulation 
in terms of uj,hu,hi, and consider, for all uj P, 

n%uj), = nu Rf{co) U Rlico), (2.1) 

where U = Mx]0, 1[, 

Rf{uj) = {x = (xi,X2), > X2 > -ehioT^^/^{uj)} , 

Pli^) = [X = (X1,X2), < X2 - I < shu OT^^/^{uj)j . 
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Figure 1: The rough domain fl^. 

We also define Eq = {0}, Si = M x {1} the horizontal boundaries of 0. 

Note that our modeling is derived from stochastic homogeneization, where domains with 
small holes are described with such random fields (see again ^S], [7j among 

others). However, we emphasize that the classical tools of homogeneization (such as compen- 
sated compactness, two-scale convergence) do not apply to our boundary problem: broadly, 
boundary layers are not seen in weak convergence processes, so that they do not allow to 
recover precise energy estimates. To prove the theorems of the next sections will require a 
precise construction of approximate solutions. 

Note also that this random framework includes the periodic one. Take P = T the unit 
torus, V its borelian u-algebra, and tt the Lebesgue measure. For any periodic function F on 
T, any yi,uj in T the formula 

fiyi,u;) = Fiyi+u;) = FiTy,iu;)) 

allows to "randomize" the periodic structure. 

In addition to the rough domain fi"^, we need to define boundary layer domains, that will 
be useful in the study of Navier's wall law. Namely, 

7^/(w) = {y = (yi,y2), 2/2 > hioTy,{Lj)} , 
T^u{uj) = {y = (yi,y2), y2 < Ko Ty^{uj)} . 

We will denote 

Hi^uiuj) = n+juj) u So u ni-joo), 

with 

7^±Ja;) = 7^^,,(a;) n {±2/2 > 0}. 

2.2 Main results 
2.2.1 Dirichlet wall law 

The first step in the analysis of system (jl.lf) , H1.2|) , is to prove the existence and unique- 
ness of solutions. For given parameters £,uj, this question has been adressed by Ladyzenskaja 
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and Solonnikov in article jl9j . We also refer to for good overview on channel flow prob- 
lems. The results of |19j yield existence and uniqueness of a solution •) in the space 



B2{n^io)) = {ue Hli^^in^iio)), sup i [ 

[ R>1 Jq 



\Vu\^ < +00 > , 



where 



n%uj,R) = n%uj)n{o <\xi\ < R} , r>o, 



for small enough flux (j) < 0o- However, we can not apply this result, as the dependence of 4>q 
with respect to e and to, as well as measurability properties with respect to lu are not clear. 
Hence, we show directly the following existence and uniqueness result: 

Theorem 1 There exists (j)Q > and eo > such that: for all cj) < (j)o, for all e < Eq, for all 
LO ^ P, system (jl.lj) . ()1.2|) . H1.3|) has a unique solution 

Moreover, if we denote u^{uj, •) the zero- extension ofu^{uJ, •) outside il.^{u;), then the mapping 
is measurable. 

As will be clear from the proof, the solution is a perturbation of the following flow: 

{u^{u!, x) = (6</)X2(l — 3:2), 0) , X € Q, 
n°(u;,x) = 0, X £ - n. 



(2.4) 



Note that in restriction to Q, is simply the Poiseuille flow, that is the solution of H1.4() . 
p.5j) . We show the following estimates 

Theorem 2 There exists C > such that, for all lo € P: 



sup 

_R>1 -K Jn^(w,R) 



sup- 

_R>1 -K JSo(i?)USi{_R) 
1 



(2.5) 



sup- 
I. R>i -K Jn(R) 



where 



Vl{r) = 1^ n {0 < |xi| < K\ , So,i(-R) = So,i n {0 < |xi| < i?} . 



Theorem 12.51 expresses that the Dirichlet wall law yields a 0{e) quadratic approximation of 
the real solution. The obtention of the L^-estimate relies on a duality argument. 
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2.2.2 Navier's wall law 

The Dirichlet wall law, that is the approximation of by does not account for the 
behaviour of near the rough bondaries. To derive a more accurate wall law, we carry a 
boundary layer analysis. We show that, for small e, is close to 

e / ^ 0/ \ , TT f , TT f xi X2-1 

Kppi^^x) =-« (w,x) +eUi luj,—, — ] +eUu [ to,— 



■ eu^{uj, x), 



e e J \ e e ) (2.6) 



where \Ju,i is an appropriate boundary layer term, and an appropriate Poiseuille type flow. 
Precisely, Vi^u = U^ui'^-.y) satisfies a Stokes system with jump conditions: 

- v/^Ui^u{^, ■) + vPi,u{^, = 0, ye nf^^iuj) 

div C/i,n(c^, •) = y ^ T^iA^), 

[Ui,u{^r)]ko=0, (2.7) 
[d2Ui^u{^^, •) - Pi,„(w, •)e2] Iso = -Gt/*, 
, Ui^u = 0, ye dni^u{(^). 

In the case of periodic roughness, the solvability of system (|2.7|) is direct. Morever, using the 
Fourier transform in tangential variables, it can be seen easily that Ui^u converges exponentially 
fast to a constant as 2/2 goes to infinity. The random setting requires much more work. We 
prove in section [l] the following 



Theorem 3 For all </> > 0, there exits a unique variational solution Ui^u of (|2.7|) . in the 
sense of paragraph \4 ■ It is almost surely a classical solution, and satisfies 



sup^Ef / |VC/i,,(-,y)|'dy ) < +00, 

R>i R \Jni,^{-,R) ) 

sup— / I VC// u(i-i;, y)|^ < 00 almost surely , 

R>i R Jn,„(uj,R) 



(2. 



where as usual 



7^^,«(w, R) = 7^^,„(^J) n {0 < |yi| < R}. 
Moreover, there exists a measurable map: Ul^^: P^R^, with Uj^ ^ = ^' such that 

sup E I 1 / \Ui^u{lo, yi,y2)dyi - [/^^(u;)^ dy 

R>1 \^J\vt\<R 



yi\<R J J'^-oo (2.9) 

\Ui,u{'^,yi,y2) - Uf°^{uj)\ > 0, locally uniformly in yi. 

The convergence of Ui^u towards a constant is proved using the ergodic theorem. Contrary 
to the periodic case, we are not able to precise the speed of convergence (see section^ for all 
details). 

We then establish the estimates on — n!:„„: 
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Theorem 4 Let Ui^, 

satisfies, as e — > 0; 



as m 



theorem\^ given in (|4.2() . The approximation u^^pp of (|2.6|) 



sup — E 
R>1 R 



sup — E 

i?,>i R 



ln{R) 



|n^(w,-) 



"app 



(2.10) 



In the periodic setting, for which the boundary layer profiles Ui^u — decay exponentially 
with 2/2) the bound o(e^) turns to O(e^) for the estimate, and O(e^) for the estimate. 
But in the general random setting, we are not able to improve the bound o{e'^): it is related 
to the speed of convergence of Ui^u as 2/2 — *■ +00. Theorem ^ extends the results of Jager 
and Mikelic, |151 theorem 1, pll3] for periodic roughness. As they consider domains with 
lateral boundaries, they need to add a lateral boundary layer term, which yields a less precise 
estimate: the O(e^) and O(e^) bounds are replaced by O(e^) and O(e^) bounds respectively. 
The unbounded setting allows to avoid this loss of accuracy. 

We may now justify the relevance of Navier's friction law. Let w^(u;, •) the solution of 

V - Vv + Vq - I'Av = 0, X G 0, 
div V = 0, X £ Q, 



(2.11) 



Vl 
cr{xi) 



4>, V2{id,-)\T:o,i = 0, 



with Navier's condition 
vliuj,-) 



+£ ai{uj, 



dvl{uj, ■ 
dx2 

dvl{uj, ■ 
dx2 



X2 = 0, ai{uj) 



X2 = 1, au{uj) 



TTOO 



(2.12) 



with f/j°°, as in theorem IHI Note that by linearity of (|2.7j) . constants ai^u{^) depend only 
on the roughness, not on the flux (p. We state: 

Theorem 5 



sup E 

R>1 H 



in(R) 



\u'{u,-)-v'{uj,-)Y 



0. 



Thus, the Navier's friction law leads to a (slightly) better approximation than Dirichlet's 
law. Again, the proof of theorem El involves the ergodic theorem, which prevents quantitative 
bounds (the o(e) can not a priori be precised). We end this presentation of the results with 
two remarks: 



1. All the estimates of section [2.2.21 involve the expectation of the spatial quadratic norms. 
Due to a lack of deterministic control of the boundary layers, we are not able to obtain 
almost sure estimates in i?2(^^)) like in theorem |21 However, one can deduce informa- 
tions with high probability. For instance, theorem El and Tchebitchev inequality show 
that: for ah R, S > 0, 



P 



n{R) 



\u%u;,-)-v'iiJ,-)\ > 



0. 
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2. Under an assumption of ergodicity on the group {t^), that is 

e V, {Th{A) = Ayh) {tt{A) = or 7:{A) = 1) , 

the constants ai^u that appear in (|2.12j) are independent of uj, as a consequence of the 
ergodic theorem (see section 0}. In such case, one does not need to know the shape of 
the boundary to determine the appropriate coefficient in Navier's wah law. It can be 
deduced almost surely from a numerical computation involving another boundary. 



3 Justification of Dirichlet wall law 

3.1 Well-posedness 

This section is devoted to the proof of theorem^ The solution {u^,p^) is searched as a 
perturbation of the Poiseuille type flow {u^,p^), where is defined by (|2.4j) . and 

p^{uj,x) = -12i^(j)Xi, Vx e ri^(w). 

We denote = + w, p'^ = p^ + q, and consider the following system (dependence on u is 
omitted to lighten notations): 

' w-Vw + u^ -Vw + w- Vit° + V(? - i^Aw = f, X (So U Si) , 

div w = 0, X e 0"^, 

MIsouEi = 0, [92-;^ - qe2] Isousi = -^(p (3-1) 



Wl 
cr{xi) 



0, w\en^ = 0, 



where 



/^(w, x) = 0, X in Q, 

f{uj,x) = {-12iy(l),0), xmn'{uj)\n. 

The proof of theorem^ is divided in three steps: 

1. We consider the linear problem 

' {v + u'^)-\/w + w Vii" + Vq- v/\w = / + div G, 
X G J]*" \ (So U Si) , div t(; = 0, x G 0% 

MIsoUSi = 0, [d2W - qe2] IsqUEi = 4> 

0, 



(3.2) 



Wl = 0, w\dne 

a{xi) 



where for all uj, 



G{u;,.) G LLc(f^'(^)) = {g, f^J ^^^.^^^ \G\^ < +oo}, 

T<\xi\<T + l 



fiu;,.)GLli,,{n'{Lo)), f = Om Q, 



(uJ, ■) G L: 



uloc\ 
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We show the existence and uniqueness of a solution 

2. Thanks to this hnear analysis, we show that, for e < Eq, (f) < (j)Q small enough, there 
exists a unique solution w of (|1.H) , thus of a solution of , as in theorem ^ 

3. We show measurability properties of . 



3.1.1 Linear problem 

To build a solution of ()3.2() . we consider the following approximate problem: 

{v + u^) ■ Vwn + Wn ■ Vu*^ + Vg„ — uAwn = f + div G, 
rr G \ (SoCn) U , div zi;„ = 0, x £ Q"^, 

[w^n]|so(n)UEi(n) = 0, [d2Wn - ^62] |i;o{n)USi(n) = 

I, Wn\dn^„ = 0, 



(3.3) 



where is a short-hand for the bounded domain Q^{uj,n). A formal energy estimate yields 



u I \VWn\ 



fw„ 



GVWn 



iWr, ■ VvP)W. 



+ Wn-fp I Wr, 

JSo{n) JEi(n) 



uloc 



IT' \\Wn\\L'2{T,o{n)UT,i{n)) 



(3.4) 



We remind the Poincare inequality 
valid in a domain 

Oa = {x, a < xi < (3, 7;(xi) < X2 < ^u{xi)}, sup {'ji - 7„| < a, (3.5) 
for all (f in H^{Oa), with 99 = on one of the boundaries 7^,1* • 

It is easily deduced from (|3.4j) 

that: for all cp G H\Q'^ \ 0(n)), 99 = on aj]^ \ 0(n), 



\\f\\L2(ni^~n{n)) < C'^^ l|V93||L2(Q|,_f^(„)), 

IIV'llL2(So{n)USi{n)) < C e'^^'^ \\Vip\\L2(_Qe^\Q|^n))■ 
We infer from ((231), that 

|V«;„|2 <c(V^(e||/||i2 +||G||i2 +^/F||<A||i2 )||Vw;„||z.2(n.) 

\ \ uloc ii.l.nr n.l.nr. I • "■ ' 

+ '/'l|Vu'„|||2(f^,) 



(3.6) 
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Hence, for (p < cpo, e < Eq with (po, eq small enough, we get (uniformly in the variable lu G P) 



1 



-^W^wnh^n^j < C e||/||i2 + ||G||^2 +^\\cl)\\^2 

\f Xt \^ uloc uloc uLoc 



(3.7) 



By Lax-Milgram lemma, this estimate yields existence and uniqueness of a solution Wn of 

(ESI). 

We now wish to let n go to infinity. However, estimate 1)3 .Tf) degenerates in this limit. To 
obtain compactness, we will follow ideas of ^Hl and localize estimate (|3.7j) in a band of width 
rj. To do so, we multiply the first equation of 1)3. 3|) by tu^, and integrate over 

VtniR) = f^^n{o < \xi\ < R}. 

Proceeding as above, we get, for (j) < cpo small enough, 



iy\VwJ^ <Cie^\\f\\l, +\\G\\l2 + eUWl^ 

uloc uloc uloc 



(R+l) 



+ 



+ 



I'diWn ■ Wr. 



x-i\=R 



+ 



Vi\U\ 



xi\=R 



+ 



Gei ■ u 



\xi\=R 



We then integrate with respect to i?, fro r/ to r/ + 1. We deduce 



F{n) <C{e 



+ \\G\\l, +e 

ULOC uloc 



2 

uloc 



+ 



+ 



Vl\Wn\ 



+ 



(7?+l 
qnWn.l 

Gei ■ Wn 



where 



/ u \Vwn\'^dR, nf^{r],r] + l) =nf,n{7] <\xi\ <r] + l}. 

By Poincare inequality (|3.4j) . and Sobolev inequality 

one has easily that 



(3. 



diWn ■ Wr. 



f1^(»7,»?+l) 



< G I \VWn\\ 



Vl\Wn\ 



Q^(»7,»7+l) 



Gei • Wr 



Ofi(r),»7+l) 



uloc^ n) 



|Vw„| , 



n|i(r,,r,+l) 



1/2 



I VWn 
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The treatment of the integral involving the pressure is exactly the same as in jl9| . Boundary 
conditions and incompressibility of Wn imply that Jj^^i^/j ^^n,! = for all R, which yields 
Jne ^+1) uJn^i = 0. It is then well-known (see p[Sl for detailed description) that: there exists 

99,, G 7^0(0^(7?, ry + 1)), with div (^^ = wi, 

and 

\\'Pv\\H^{n;i,{r,,v+l)) < C'lkl||L2(nfi(r,,r?+l)), 

where C is a positive constant independent on e, n and ??. We can write 

qnWn,l = / gndiv ip^ = - Vg„ • ifr,- 

Using the expression of Vg„ in and after several integrations by parts, 



qnWriA 



This leads to the inequality 



< 


In 







qnWn,l 



<C l + ||u| 





+ 


j 








Jn 








+ 


[ (/ + div G)Wn 









uLoc 



\VWr, 



+ C{e^\\f\\l, +\\G\\l, 



Together with previous bounds, this yields 



|Vtt;„|| 

+ C(e''\\f\\l, +\\G\\l, 

uloc uloc 



ULOC 



that is 
where 



F{ri)<CiF'{i-i)+ C2(l + ??), 



(3.9) 



(3.10) 



Ci = Ci 1 + \\v 



uloc ' 



+ \\G\\'i, +d'^"2 



C2 — C2 -r in-T||^2 -rt||(/^||^2 



This last equation is a reverse Gronwall type inequality. By estimate ()3.7() . up to take a larger 
C2, we can suppose that 

F{n) < C2{Ci + n + l). 
From (|3.1()j) . we deduce easily that for all r/, 

+ 



e ci "i?(n) - e '^i '^F{r]) > -C2 



Ci 



e ci dt 



>C2{Ci+n + l)e Ci"-C2(Ci + r? + l)e Ci'?. 
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It follows that for all rj, 



F{^)<C{s'\\f\\l, +\\G\\l, +e\\4>\\l, )(ry + l) 

^.7„„ uloc 



and finally 

supi/ |Vu;„|2<c(e2||/||2, + ), (3.11) 

77>1 ''7 J Q.^ {t}) uloc uloc uloc/ 

where C can be taken affine in H^jHj;^! . 

uloc 

On the basis of such estimate, one can, for all w, extract a subsequence w^^f^^) that 
converges weakly in H^^^ {^"^ (u))) . The limit w{io,-) is a solution of (|3.2j) . and satisfies the 
estimate 

supi/ |V«;p<c(6^||/||i. +||G||i. +e||<^||i. ). (3.12) 

^>]^ ?y (^r)) uloc uloc uloc/ 

It remains to show that w is in H^^^^. The argument is almost the same as in |19| p745]. 
Let r > 3/2. By energy estimates carried on the translated domain 

n''{T,R) = n'{R) + (r,0), 

we obtain similarly to (|3.1fl|) 

Friv) < CiF^iv) + C2(l + V), Friv) = [ \Vw\^ dR. (3.13) 

Jri Jq~{t,R) 

By (|XT^ . we have 

FAT)<c(e'\\f\\l, +\\G\\l, +eml, )(r + l), 

\ uloc uloc uloc/ 

and reasoning as above, we get: for all < r/ < r, 

F.(r?)<c(e2||/||2 +||G||2 + ,||^||2 \ + i). 

\ uloc uloc uloc/ 

Applying this inequality with r] = 3/2 yields 

ll^llk <c(e'\\f\\l, + \\G\\l, +eml. ) (3.14) 



For uniqueness, one must show that the solution w of 1)3. 2() with f = 0, G = 0, (f) = is 
identically zero. Inequality (|3.10|) turns into 

F{7]) < CiF'{rj), Wi] > 0, (3.15) 

from which it follows easily that 

limsupF(r7)e~'^i'' > F(0). 

r;— >+oo 



Together with (|3.12|) . we deduce w = 0. Note that uniqueness does not only hold in 
1 ( 

uloc\ 



H^i^^{Q^ {uj)) but in the wider space B'^in^uj)), for any uj £ P. 
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3.1.2 Nonlinear problem 

This section is devoted to the well-posedness of (|;il|) . for any fixed lo & P. Therefore, we 
consider linear equations (|3.2() with special choice 

f = G = 0, and = -6(p. 

Note that ||/'^||j;^2 < C^/e. We define the application 

uloc 

A : HiU^'iu;)) ^ hIU^'{u;)), v{u;, •) ^ w{o,, •), 

with w{lo,-) the solution of ^'A.2\i . Then, for e < eq small enough, A is a contraction in 
restriction to the unit ball of Hl^i^^{Vf{uj)), . 

Indeed, let , in the unit ball of //^^Q^(ri^(u;)), and denote = A(u^), w'^ = K{v'^). 
By estimate 1)3. 14() . they are bounded by 

1^1' OC \ ULOC uloc / 

for e < eo small enough. The difference w = w"^ — vlP' satisfies (|l-{.2j) . with 

v = v^, f = Q, G = {v'^ -v^)®w'^,^ = 0. 
Again, using (|3.14j) leads to 



kllm <C||(t;2-t,i)0u;2||^ <C\\v^-v^\\hi W\\h\ 

uloc ULOC ULOC ULOC 



(3.16) 



1 

uloc 



Hence, for e < eo small enough, A is a contraction in restriction to the unit ball. By the Banach 
fixed point theorem, we deduce the existence of a unique solution w{lo, •) of 1)3. 1(1 in the unit 
ball of Hli^^{Vt^{uj)), for all uj £ P. Back to the original variables, there exists a solution 
u'^{uj, •) of p.lj) . which is unique in the ball of center and radius one in i^^j^^(Q^(a;)). 

Uniqueness in the space i?2(^^('^)) is deduced easily from T^, theorem 2.3, p739]. The 
idea is still to obtain local estimates on the difference of two solutions. Due to the quadratic 
term in Navier-Stokes equations, the inequality 1)3. 15(1 is modified by a nonlinear term, but 
still leads to uniqueness. We do not give further details, and refer to ^Hl- 

3.1.3 Measurability 

To conclude the proof of theorem it remains to check the measurability properties of 
a; 1-^ n^(a;, •), P i— > Hj^^i^), where the ~ stands for the extension by zero outside rj*". We 
know from the previous section that = tt^ — is the fixed point of a contraction. Precisely, 
for all u; G P, w{uj,-) is the strong limit in H^i^J^^¥{ui)) of the sequence {w^{ui,-))kizf^, 
satisfying 

xGrj^\(SoUSi), div = 0, xeQ', 

[w''+^]koU^, = 0, [52^; - qe2\ |eoUEi = -60 (3-17) 
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in which we take llu^'^ll/^i < 1, with uj w^iuj^ ■) measurable. Thus, it is enough to show 

uloc 

that for all k, u> ^ w^{uj, •) is measurable from P to Hi^^{M.'^). Therefore, we will prove that 
for any function 

v{u>, •) G H^ig^{Q'^ {uj)) , u! I— > v{uj, •) measurable, 

the solution w of (|3.2j) (with f = , G = 0, cj) = —Qcp) is such that uj i— > w{u), •) is measurable. 

With previous arguments and notations, for all w, there is a subsequence 'w^^(n)i^r) of 
Wn{uJ,-) that converges weakly in Hj^^iuj'^ {oj)) to w{ijJ,-). By uniqueness of w as a solution 
of (|.S.3j) . one can easily check that the whole sequel Wn{u), •) converges to w{u), •). Hence, by 
Pettis theorem it is enough to show the measurability of a; i— > Wn{uJ, •). 

To do so, we follow ideas of [Ij on an elliptic problem from homogeneization. Let us 
introduce 

Cn=]-n,n[x]-2,2[, 

and the spaces 

Vn = {we Hi {Cn) , div u; = 0} , 
Vn{oj) = {w eVn, w = Q\nCn \ ^"{uJ, n)] . 

Let vrn(w) : Vn ^ Vn the orthogonal projection on Vn{oj)- By definition of {P,V), it is easy 
to show that for all n, the set-valued map 

J'n-. {P,r)^Vc{H\R^)), c^^KM, (3.18) 

is measurable. Following ^01 , we remind that a set-valued map 

from a measurable space P to the non-empty complete subsets of a separable metric space X 
is measurable if: for all open subset O of X, the set 

{w, J^(aj) n O / 0} 

belongs to V (in our case, this set is open, and so belongs to the u— algebra V). For all w 
in Vn, there exists by property (|3.18|) a sequence of measurable selections crj{u)) such that 
(Tj(a') — > Trn{io)w strongly in Vn {c.f. \1(A theorem III. 9, p67]). We thus get, for all w, w' in 

Vn, 

{{'Kn{uj)u,v)) = \\m.{{aj{Lo),w')), 

3 

where ((•, •)) denotes the usual scalar product on Hl{Cn). Thus, for all w, w', 

w 1-^ ((7r„(a;)t(J, w')) 

is measurable, which by Pettis theorem yields the measurability of 7r„. 

Assume now that uj v{uj, •) is measurable. There exists q„ > 0, such that for all lo, the 
solution Wn{LO,-) of H3.2|l is the fixed point of the contraction: 

F{uj) -.Vn^Vn, vr„(u;)( an(A{u)w - l{uj)] +w], 
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where ^(a;) G Vn, and A(lo) : Vn Vn are defined by 

(( A{lo)w, w' )) =v / {Vw ■ Vw' + {v{uj, •) + n°(w, •)) ■ Vw ■ w' 

J Cn 

+ w ■ Vn°(u;, •) • w'^ 

Thus, Wn{uJ, •) is obtained as the limit of the sequence 

<+i = F(a;)«). 

It is clear that uj ^ l{uj) and lo ^ A{u)) are measurable. Together, with the measurability 
of 7r„(a;), it follows easily that lo ^ F{uj)w is measurable for all w in H^(M.'^), hence the 
application lo ^ Wn{LO, •), from P to H^(M?), for all j. As m„ is the limit of Wn, theorem^is 
proved. 



3.2 Estimates for Dirichlet wall law 

We now turn to the proof of theorem [21 Instead of estimates in B2 norm, we will show the 
following more precise L^i^^ estimates, valid for (j) small enough: 

m<CV~e (3.19) 

ll-II^L.(^oa) ^ (3-20) 
in)<Ce. (3.21) 

In fact, the first of these inequalities is just (|3.14|) with G = 0, and the second one is an 
easy consequence of (|3.6() . So we just need to focus on the third inequality, which will 
be proved through a duality argument. The main task is to estimate ||w||j;^2(f2^) (where 
J7r = n {— i? < xi < i?}). In fi, we can see w as the solution of a modified Stokes system, 
with inhomogeneous Dirichlet conditions on Sq and Si : 

' Vg - vl\w = -{w + u^) ■ Vw - w ■ Vn°, x G 

< div u; = 0, xefl (3.22) 
^ w given on So.i, 

the boundary values of w being controlled by (|3.2()j) . We introduce the adjoint equations in 
i^R, that is the Stokes system with vanishing boundary conditions and non-zero source term: 

' Vvr — I'Av = f, X G 

< div u = 0, xe Qr, (3.23) 

, v\dnR = 0, 

where (p £ L'^{0,r). We proceed in two steps. First, we establish two regularity lemmas for 
(|3.23j) . Then, by appropriate choice of 93, we obtain the claimed estimate (|3.21|) . 
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3.2.1 Two lemmas for the Stokes system 

It is well known from the elliptic regularity theory that the solution {v, vr) of (|l'{.2i-{j) lies in 
X as soon as 99 G L^, and that there exists C(Or) such that ||f ||//2(q^) + < 
C(J7/j)||(/?||j;^2(Q^) (here vr is the zero-mean determination of the pressure). However, the 
constant C{^}ji) depends a priori on R. We need an estimate independent of R, which may 
be obtained by localization. 

Lemma 6 There exists C > such that for all R> 1, the solution (f,vr) of (|.S.23|) satisfies 

hWH^injt) + l|V7r||i2(n^) < C\\(f\\L2(^Q^y 

Proof. 

Let X e C~(M) such that x = 1 on [0, 1], x = outside [-1, 2]. Let Xk{xi) = xi^i - k). 
We will also make use of two other truncation functions 9l et On such that 6l ^ C°°(M), 
6»i = 1 on ] - 00,2], 6Il = on [3,cx)[, and 6lj?(xi) = Oii-xi). Let K = [R\ - 1. We 
localize the equations thanks to the functions 9i, Xk and Or: let V].{xi,X2) = Xk{xi)v{xi,X2) 
{—K < k < K — 1), vl = Ol{xi + R)v, vr = 6r{xi — R)v. The function Vk (resp. vl, vr) 
vanishes outside the domain = {k — 1, k + 2) x (0, 1) (resp. Dl = {—R, —R + 3) x (0, 1), 
Dr = {R—3, R) X (0, 1)). Let Ck be the mean of the pressure vr on Dk- We set vr^ = XkiT^ — Ck)- 
We define similarly vr/,, 'kr. The functions Vk,TTk solve the following equations in : 

' Vvrfc - uAvk = (fk, X £ Dk 

< div Vk = ^k, X £ Dk (3.24) 

. VkldDk = 0, 

where 

^k = XkV - i^(2Vxfc -Vv + v ■ Axk) + Vxfc(vr - Cfc), 
7fc = Vxfc • V. 

In Dl and Dr, vl and vr solve analogous systems (note that vl equals zero on the left 
boundary of Dl , for v does). The classical elliptic theory applied to this system yields 

ll^fclli2{Dfe) + l|V^fcllL2(Dfc) ^ C'(ll¥'fclli2(£,^) + hklln^n^))- 
Of course, the constant C is independent of k. We have 

K-l 
k=-K 

+ C^h-Ck\\l2^o,)- 

k 

We can control vr — in the following way. One can find C > such that for all / G L'^{DjS) 
with / = 0, there exists a function ip G HQ{Dk) satisfying div ip = f and ||V'll_f/i{Dft) — 
C||/||l2(£)^). Hence we have for such / and ip, 



/ (tT - Cfe)/ 




/ Vvr • V 









<c||v^IIh-(d.)II/IIl^(d,). 
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Thus we obtain 

h - Cfc||L2(Dfc) < C'l|Vvr||^-i(£,^) < C{\\ip\\L2^Dk) + l|Vv|lL2{Dfe))- 

We have similar results for vl and vr. Noticing that on (k, A; + 1) x (0, 1), we have = v 
and VvTfc = Vp, and similarly with vl, vr on {—R, —R + 2) x (0, 1) and {R — 2,R) x (0, 1), 
the previous inequalities lead to 

K 

k=-K-l 



(where indices k = —K — 1 and k = K stand for k = L and k = K to lighten notations). 
Last bound comes from Poincare inequality and from the standard energy estimate of H3.23() 
in the whole domain This ends the proof of the lemma. 



Our second lemma is an H^i^^ estimate for the Stokes system in Qr, uniform in R. We 

set 

||f ||?2 ^ = sup / \vf 

uloc\ Rl -R<t<R-iJ{t,t+1)x{0,1) 

and ||f \ = ll^^ll?2 \ + l|Vf ||r2 \- 

II ''H'^, (Qr) II IIL^, (Qr) II ll-L , (^r) 

Lemma 7 There exists a constant C > such that for all R > 1, we have 
Proof. 

The proof is very similar to the calculations which led to the estimate (|;-i.l4|) , so we only 
sketch it briefly. We write an energy estimate on the domain {—R,r) x (0, 1), and integrate 
with respect to r, from 77 to r/ + 1, to obtain a backward Gronwall type inequality: 

F{7j)<c[ \ip\^ + CF'{r^)<CM\l, {R + r^ + l) + CF'{rj) 

J{~R,-n+i)x{o,i) "'""^ 

{—R < rj < R — 1), where 



2 ! 

uloc 



F{r]) = / \Vv\'^ dxdr. 

Jti J{-R,r)x(0,l) 

For rj = R — 1, we have 

FiR - 1) < ||V^;||i2(^^) < Mht^n,) < 2R\Ml 
thus we obtain 

F{rj)<CM\l, (R + V + l)- 

uloc 

Similarly, we have, for —R<rj<R — l, 

G{r])= f \Vv\^ dx dr < C\\ip\\l2 (R-v)- 

Jri J{r,R)x{0,l) "'"'^ 
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Now let T G [—R + I, R — 1]. Again, an energy estimate on (r — r, r + r) x (0, 1) leads to 



\Vv\^ dx dr < C\\ip\\l2 {rj + 1) + C H'^{r]) 

r] V{T-r,T+r)x (0,1) "'"'^ 



Hr{ri) 

(0 < r/ < i? — 1 — |t|). If, for instance, r > 0, we have 



H^{R-1-t)< I 



{2t-R,R)x{0,1) 

< G{2t-R-1) 
<CMl2 (2i?-2r + l), 

uloc 

hence we conclude that for allO<r/<i?— 1 — r, 

Hrir])<C\Ml, iv + l). 

uloc 

If r < 0, we obtain a similar conclusion using F instead of G. For t] = 1/2, we get 



(T-ir+i)x(0,l) 



iVwp < C\\ip\ 



uloc 



The result follows. 



3.2.2 The duality estimate 

Now we estimate w using a duality method. Let R > 0, and let v^, vr^ denote the solution of 
the Stokes problem 1)3. 23(1 in 0,^ with source term if = Using these auxiliary functions, 

let us show the following inequality for (p small enough: 



< CRe^ + C 



\w\ dx2. 



(3.25) 



xi\=R 



To prove this, our starting point will be: 



w 



w ■ (-uAv^ + Vtt 



/ {-uAw -v^ + w- Vvr^) - / 



uw 



QyR 

dn 



(3.26) 



We have to control the various terms appearing in (|3.26j) . We recall the equation satisfied by 
w in 17 ij: 

-vAw + {w ■V){w + u°) + u\diw + Vg = 0. 
Thus integrations by parts lead to 



vAw ■ 



w- {w- V)v^ + -{w- V)u" - u^w ■ div 



(3.27) 
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Let (ao, un) be a regular subdivision of the interval [—R, R], such that oq = —R, aat = R, 
and 1 < afc_|_i — flfc < 2 (hence < 2i?). We divide the domain into the corresponding 
cells Ck = (afc_i,afc) x (0,1). Then we can write, using the Gagliardo-Nirenberg inequality 
in each cell Ck, 

N 



LHCk) 



Wv 



R\ 



lL2(Cfe) 



<Y.C\\w\\l2^c,)\MhHCu)\\'^A\lI,^^ 

k 

<Y.^^\\M\lHc,) + C{A^r'\\w\\\, \\vXh^ ) 

' ^ ' uloc uloc 

k 



n^p) ^ ' (3.28) 

^ uloc 

thanks to the second lemma, with r] a small parameter to be chosen later. We obtain, using 
the energy inequality for ()3.23() . estimate p.l9() and the previous inequalities ()3.27|) . H3.28() : 



uAw ■ 



<iv + Ccl>)\\w\\l,^^^^ + CRr^' 



(3.29) 



Next, we treat the boundary integral in 1)3. 26() . We have, for i = 0, 1, thanks to 1)3. 20() and 
lemma 



vw ■ d2V^ dxi 



MR) 



< 



l^lkllL2{S,(iJ))l|Vt'^||L2(s^(^)) 



< CeVR\\v^\\H2 



< Cri-'^e'^R + ri\\w\\l2 



(3.30) 



We also write, using again a trace theorem, 

,c( 



vw ■ div^ dx2 



\xi\=R 



\w\ dX2 



xi\=R 



<Cv ^ \w\^ dx2 + v\\w\\l2/n„) 

J\xi\=R 



(3.31) 



(here the integrals on the sections xi = ±R are restricted to the portion of each section 
contained in O). It remains to deal with the integral involving the pressure in (|3.26() . We 
use again the subdivision of Jl/j into the cells Ck- Let Ck denote the mean of tt^ on Ck- We 
integrate by parts w ■ Vtt^ = w ■ V{tt^ — Ck) on each cell Ck- 



[ ^„ . vvr^ = V / 

= E 



Ein9Cfc 



W2{tt^ - Ck) dxi 



w ■ n{-K — Ck) 

■W2{tt^' - Ck) dxi 

- / wi{tt^ - Cl) dX2 + / Wl{-K^ - Cn) dX2 

Jxi=—R J xi=R 

n— 1 „ 

+ ^/ Wi{ck+1 - Ck)dX2- 

k=l ■Jxi=af, 



(3.32) 
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We bound the quantities appearing in the r.h.s. of (|3.32|) . We have, for i = 0, 1, 



k 

<^(C7?-V + r?||Vvr«||i2(c,,)) 



using (|.S.20j) . the Poincare-Wirtinger inequahty apphed to the zero- mean function vr^ — Ck in 
Cfe, and lemma El Similarly, we have for the two terms corresponding to xi = zizR: 



w^{tt^ — Ci^n) dx2 



xi=±R 



<C'n \w\'^ dx2 + v\\w\\h(nR)■ 

J\xl\=R ^ ' 



To deal with the last term of ()3.32|) . we use the fact that Xt((jj.) dx2 = 0, where (y{ak) 
{(afc,X2) € (here we take into account the exterior Q!^\Q oi fi). Hence 



1/2 



uloc ^ ' 

(we can apply the trace theorem in OF n {a^ < x\ < Cfc+i} for instance, with a constant 
independent of k since the boundaries of 0^ are X-Lipschitz continuous). Moreover, we 
obviously have |cfe+i — Cfc| < C'||V7r^||^2((7^y(7^_|_^), thus estimate ()3.19() and lemma El again 
imply 



/ W\ dX2 




1 Wi dX2 


Jxi=ak 








<-[ 


e / 









E 



'W^l(Cfc+l - Ck) dX2 



xi=ak 



< J^<^e||V7r-^||L2(t7,uCfe+i) 

l|2 



< C'q'-^e^R + r]\\w\\\2 



Up to this point, the r.h.s. of H3.32() can be controlled as follows: 

! wVtt^ <Cri-^e'^R + 4:ri\\w\\l2(^ ) + C [ \w\'^ dx2. 

JUr J\xi\=R 

Then inequalities (jSHH), (HOni) . (nOB and enable us to write, from (jH^Sl): 



(3.33) 



/ \w\^<{Ccl> + 8r^)\\w\\l2^^^) + CRv-'e^ + Cii-' [ 



\w\'^ dx2- 



We take t] = 1/10 and (/)< (l)o such that C0o < 1/10. This leads to 

Let z{R) = \\w\\'L2(^Q^y Inequality (pT^ can be rewritten z{R) < C Re"^ + C z' (R) . Follow- 
ing ^ni; we deduce the estimate p.5j) from this differential inequality and from the fact that 
z{t) has sub exponential growth at infinity (Poincare inequality and (|3.19|) imply that it has 
linear growth). Indeed, let z{R) = CRe^ +C'^e^; we have z'{R) - z'{R) > ilM^lM. Hence if 
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there existed Rq such that z{Ro) - z{Rq) = yl > 0, we would have z{R) - z{R) > Ae'^^'^^^/^ 
for all R > Rq, and z would grow exponentially fast at infinity, which is not the case. Thus 
we obtain for ah R> 1, z{R) < CRe^ +C^£^, that is (j^ . 

Finally, we show 1)3.21(1 . It is a direct consequence of what we have already proved. For 
i? = 1, we have ||^^||/,2(f2^) = -2(1) < Ce^ . Now we can repeat the whole proof but with another 
origin: take xi = r instead of xi = as origin, replace with Q.{t,R) = + (r, 0). This 
yields Zt-{R) = ll''^lli2(Q(^ /jj) ^ CRe^ +C^e^, with the same constant C as before. Hence we 
have again ||^w||/,2(q(,- < Ce^, and we conclude that ||w||^2 < Ce^. 

4 The boundary layer analysis 
4.1 Formal expansion 

To obtain a refined wall law at di}, one must clarify the structure of the flow near the rough 
boundary. Namely, we will show that behaves as 

Xi X2\ / Xi X2 - 1 



uj, — , — j + eUu iuj, 

+ eu^{uj, x). 



£ £ J (4.1) 



In this formal expansion, Ui^ui^, u) are boundary layer profiles, expressing the strong gradients 
of near the boundary. They are defined on the rescaled domains 7^/^u(w). They should be 
localized near the boundary, and cancel the jump of the normal derivative of u^. As usual, to 
derive the system they satisfy, we plug expansions ()4.1() in equations (|1.1|) . and collect terms 
of order e~^. Formally, this leads to the Stokes system ()2.7|) a.b. Jump conditions ()2.7|) c.d 
ensure the regularity of the approximation through dQ. Finally, (j2.7|) d is in agreement with 
the Dirichlet condition (|1.2() . 

As mentioned in theorem 13 we will show that Ui^u{<^i ■) converge to some constant fields 

as 2/2 +00. Because of these fields, which do not match the required flux and boundary 
conditions, there is an additional corrector n^(a;,x). It is given by 

' u\u:,x) = {3{U^, + U^,)xl - {m^, + 2[/~i)x2 - o) , x G 17, 
< uHuj,x) = -U^, X2<0, (4-2) 
^u\u;,x) = -Ur, X2>1. 

Note that in the interior domain Q, satisfies 

u° ■ Vu^ + ■ Vii" - uAu^ + Vp^ = 0, X G ft, 
div = 0, X € 0, 

n%, = -U^ n^ls, = -[/r, [ n\ = -{Ul^,+U^,). 
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4.2 Well-posedness of (|T7jl 

We now justify the previous formal computations, and prove theorem IHI We first consider 
the weh-posedness of system ()2.7() . In the study of Dirichlet wall law, the Poincare inequality 
(|3.4j) . applied to the channel Q^, allowed for deterministic reasoning. In the boundary layer 
domains lZi^u{'-^) that are unbounded in every direction, this inequality does not hold, and we 
shall use the probabilistic modeling to solve (|2.7)) . This approach borrows to homogeneization 
problems (see ^1|71IH]), but strong changes are needed to account for the anisotropy of our 
boundary layer domains. The general idea is to construct a solution of the type 

= V^,„ (rj,i(a;),y2 ± o Tj/i(u;)) , V^,n = Vi,„(w, A), (4.3) 

where Ty^ and hi^u were defined in subsection 12.11 This requires to introduce adapted func- 
tional spaces and variational formulations. As the lower and upper systems are treated in the 
same way, we will focus on the lower one, and drop the subscript "1" for brevity. 

4.2.1 Stochastic derivative, Convolution 

Let y : P X M M", n > 1, measurable. We call a realization of V an application 
:M2^M", y^V{Ty,{uj),y2 + hoTy,{uj)), oj e P. 

Notice that 

Rr,Suj)[V]{yi,y2) = Ru.[V]{yi + h,y2), yh,y. 

We say that V is smooth if, almost surely, is smooth as a function on M?. For all 

smooth functions V, we define the stochastic derivative of V as 

d^V{u;,X) = dMV] {yi = 0,^2 = A - h{uj)) 

We also introduce 

dxV{u,\) = d2R^[V] {yi = 0,2/2 = A - hiuj)) 

and the gradient V^^ = {daj,dx)- Notice that, almost surely, 

Ru. N^V] = VR^[V] (4.4) 

which allows to define for smooth V partial derivatives at any order d^V. 

For interesting examples of smooth functions, one can use a convolution process. Precisely, 
let p = p{y) G M a smooth function with compact support. For all V G L\^^ (P x R) = 
(IxL^iP X K), K compact, we define 

p*ViLO,X)= piy)V {Ty,{uj),X + y2 + hoTy,iu;) -h{io)) dy. (4.5) 

The following lemma is easy and left to the reader. 

Lemma 8 Assume 1 < p < oo and V ^ LP {P x M), resp. Lf^^. Then p*VgLP{Px M.), 
resp. L^g^, and is a smooth function. Moreover, the following identities hold: 

Ruj[p*V] = P*RLLi[y], o-.s. in IV, 
where p{y) = p{—y). Besides, ifV is smooth and V^^y G Lj^^(P x M), then 
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4.2.2 Functional spaces 

We can now introduce functional spaces adapted to (|2.7|) . Let T)q{P x R^) the set of smooth 
functions Lp = ^{oJ, A) € M^, satisfying 

i) For all a, [5, d^ip G L~ (P x M) . 

ii) Almost surely, the support of is included in {y, \ > y2 + h o Ty^{u>) > 6} for some 

positive constant 5 = 6{uj, ip). 

Note that realizations of functions of T)q{P x M_|_) have a support that is strictly included in 
TZ{uj), almost surely in lv. We clearly define a scalar product on T>o{P x ]R_|_) through 



JPxM. 

Let Dq{P X M+) the completion of 'Dq{P x M+) for the norm Hv^Hdq = ((v^)'/^))^- It is of 
course a Hilbert space. We state the following properties: 

Lemma 9 We have 

i) Do{P X M+) ^ CO(M+; L\P)) ^ Ll^{P x 

ii) For any V in Dq{P x ^/lere exists a unique D^V £ (P x 1R+) such that 

f (D^V) ip = - [ (V^ip) V, V99 G Po(^' X (4.6) 

JPxM+ JPxR+ 

Moreover, \\D^V\\ = ||V^||do- With a slight abuse of notation, we shall write V«, = 
{dx,di^) instead of D^^ ("weak stochastic derivative "). 

iii) For any V in Do{P x almost surely, Ruj[V] G HI^^{TI{uj)) , -Rw[l^]|a(7e(a;)) = 0, 

VR^[V]=R^[V^V]. 

Proof. 

i) We can argue for V £Vq, as the result for Dq follows by density. Almost surely, 
belongs to Hl^{n{uj), R^[V]\9n{cu) = 0. Thus, 



/ \V\\u;,X)dP = 11 dyi [ \V\^{Ty,{u;), X)dP 
Jp ^ -^{-1,1) JP 

= \\R^[V\[yi,X-hoTy,{u;))\\l,^_^^^^ dP 

< C{X) ||Vi?^[y](x,z)||i.(^(^^,)) dP 

= 2C{X) [ \V^V\'^dPdX 
where the inequality comes from classical trace theorem and Poincare inequality. 



25 



ii) Let V £ Dq, and Vn £ Vq, converging to V. As ||K-Ki||do = \\^u;Vn-VujVm\\, V^K 
is a Cauchy sequence in L^{P x thus converging to some D^^V with ||Dt<;^||L2 = II^IIdq- 
As 

/ (V^K) ip = - [ (Vu^if) Vn, yip G VoiP X M+), 

JPxR+ JPxR+ 

relation (|4.6j) fohows (notice that the right hand side converges by i)). 

It remains to show uniqueness of D^^V. This fohows easily from the following general 
result: for any W in {P x M+), resp. Lf^^, 1 < p < oo, there exists W„ G ^o(-P x 1^+) such 
that 

Wn > W in L'P {P X M+) , resp. Lf. 

n— »+oo 

To see this, first notice that functions of L°°{P x M+) with compact support in A are dense 
in LP{P X M+), resp. L^^^. We therefore assume that W is such a function. We then proceed 
by convolution: let pn an approximation of unity, and set Wn = Pn * W. From lemma |H1 we 
deduce easily that Wn G 'Dq{P x M_|_). To obtain the convergence, we involve as above the 
realizations of Wn and W, by extra integration with respect to yi. Namely, for W £ L^, 



I dP\Wn-Wf = l[ dyi I dP\Wn-W 

JPxR ^ -'(-1,1) JPxR 



P 



'(-1,1) 

<\l I dy\R^[pn*W][y)-R^[W]{y)f 

IP J(-l,l)xR 



2 



< f dP I f dy\pn* Ru [W] {y) - R^[W]{yW 

JP Vj(-1,1)xR , 



'(-l,l)x 

By standard results on the (standard) convolution, the integral inside the parenthesis con- 
verges to as n — > +00, almost surely. The dominated convergence theorem allows to 
conclude. Similar reasoning holds for W S L^oc^ with P x K, Kcompact, instead of P x M. 
iii) follows again from a density argument, and is left to the reader 

We finally introduce the (closed) space of divergence-free vector fields: 

VoiP X ^+) = {Ve Do{P X M+), V^-V = d^Vi + dxV2 = 0} . 

4.2.3 Variational formulation 

The keypoint to solve (|2.7|) is to search the solution U{u, •) as a realization: almost surely, 

U{u, •) = R^[Vl ye Vo{P X M+). (4.7) 

At a formal level, if we substitute (|4.7|) in (|2.7|) . and test again R^[ip],Lp G Vq{P x M+), we 
end up with 

y I V^y : V^v9 = £(99), e Vo{P x M+), (4.8) 

JPxR+ 

where the linear form 

C{lp) = j ^{uJ,+h{uj))dP 

comes from the inhomogeneous jump condition on the normal derivative. Note that C is well- 
defined and continuous on Dq{P x M+) by lemma IHI This suggests the following definition: 
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U is a variational solution of ()2.7|) if it satisfies H4.7I) . ()4.8() . 

Notice that by Riesz theorem, there exists a unique solution V to ()4.8|) . and so a unique 
variational solution ?7 to (|2.7|) . From lemma El one has almost surely, 

?7(a;,-) G ^'Wc^)), f/(w, OlaCT^H) = 0, VC/(l^, •) G (7^(L^, i?)) , Vi?, 

with for ah i? > 1, 

1^^^ |VC/(-,y)pd2/^ = llV^Flli. < oo. 
Moreover, by the ergodic theorem, we have 



R 

exists almost surely, so that 



hm ^ / \VU{;y)\^dy 
lim 4 / dy I \V^V{Ty{uj),\)\''d\ 



sup— / |Vf7(-, y)p(iy < oo almost surely . 

It remains to show that U is almost surely a classical solution. This is a consequence of 

Lemma 10 Almost surely, div U = in the weak sense, and for every -0 G {TZ{uj)) with 
div ^p = 0, we have 



V VU:Vi' = 6(j) v(yi,o)d?/. 

Proof. By lemma iii), div [/ = in the weak sense almost surely. To recover the Stokes 
equation in the weak sense, we reexpress the variational equation (|4.<S|) as 

1^ [ V^y -.V^ip + [ V^W : V^if = V(/p G Do{P x M+), 

JPxR+ JPxR+ 

where 1^ E V;j^(P x M+), which means that 



V^W : V^V = 0, yV€ Vo{P X 

PxM+ 



We consider again an approximation of unity p„ = Pn{y), and set V = pn * {'^^tv'), ^ G 
Po(^ X M+). We compute 



JPxR+ JPxR+ ^ 

= [ (Pn* (V^VF)) : V^V^^ = - [ (Vi • (V^)2(p„ * W) 



if. 
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so that • {Vujf{pn *W) = Q. Applying R ijj we deduce that 

■ /\{pn*Ru.[W]) = Q. (4.9) 
Back to (j4.8|) . we take = p„ * 99 as a test function. One has easily the following identities: 

j v^y : v^(^„ = - j {{v^f {pn * y)) • 

j ^u^W : V^(^„ = - j {{V^f {pn *W)y^ 
C{(pn) =6(t>- j (^Pn * ((t^, A) ^ (5(A - j • v?, 

where 6 stands for the Dirac measure, i.e. 

Pn* ({u!,X) '-^ S{X + h{uj))) = / dyipn{yi,-X + h{u!))dyi. 
We thus deduce from (|4.8|) that, 

-l^iV^f {pn * V^) - (V^)' {Pn *W)=Qct>-pn* ((c^, A) ^ 5{\ - h{uj))) , 

which implies that, almost surely: 

-z^A {pn * U) {uj, •) - A {pn * Ru;[W]) {OJ, ■) = 6cP ■ (pn * 6y,=o) . (4.10) 
From H4.9() . we can express 

A {pn*RAW]){uJ,-)=^P^, 

for some smooth scalar field Pi_j. Hence, if we multiply (|4.1fl|) by a divergence free vector field 
ijj G C^(TZ{uj)) and perform integration by parts, we obtain 

V(/9„ * [/) : VV = I / {Pn*tp){yi,0)dyi. 



The limit n — > +00 gives the result. The fact that U{uj.-) is regular follows from standard 
ellipticity properties of the Stokes operator. This shows the well-posedness of the boundary 
layer system (|2.7|) . 



4.3 Convergence at infinity 

This section is devoted to the last part of theorem 01 that is convergence of Ui^u as y2 +00. 
In the periodic setting, the convergence is exponential, as can be seen easily from a Fourier 
analysis. Namely, if the roughness has period T, it is straightforward that 

for some constant a independent of T. In this case, Uj^ = (^U^^^, 0^ is just the average 
of Ui^u with respect to yi. We stress that the convergence rate a/T goes to zero as the 
period of the roughness T goes to infinity. Thus, in the random setting (in which, broadly 
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speaking, all periods are involved in the roughness), this analysis falls down. We will show 
the convergence using ergodicity properties of Ui^^- For the sake of brevity, we will only treat 
the lower boundary layer Ui. Similar reasoning holds for the upper one. 

We first establish a representation formula for Ui, in terms of the double layer Stokes 
potential. 



Lemma 11 Almost surely, the solution Ui of (|2.7j) satisfies, for y2 > 0, 



where 



Proof. 



Ui{uj,y) = 
G{yi,y2) 



G{yi - y'i,y2)Ui{uj,y[,0)dy[, 
2y2 



(4.11) 



vr 



+ yiy 



Vi yiy2\ 
ym yl ) ' 



Let us denote by Ui{uj,y) the right hand side of (|4.1H) . We remind that 

Ui{uj,y) = V{Ty^{uj),y2 + hioTy^{uj)) = V{Ty^{uj),y2). 

As a consequence, 



Ui{u;,y) = V{Ty^{uj),y2), 



where 



V{^,y2) 



G(yl,y2)V(T_,,(^),0)dyl. 



Thus, Ui and Ui are both stationary with respect to yi. Moreover, they are both smooth 
solutions of the following Stokes problem in a half space: 

- AU + Vp = 0, y£ Rl, 
div [/ = 0, ye M+, 

. U\y2=0 = Ul\y2=0- 

For more on the double Stokes layer potential, we refer to ^Sl- A formal energy estimate on 
Ui — Ui yields, using stationarity: 



E 



{~R,R)xR+ 



ViUi-Ui) {■,y) 



dy 



0, MR > 0. 



One checks easily that VC/^ and Vf// have appropriate integrability, so that this estimate is 
rigorous. We deduce that Ui = Ui almost surely. 



We now introduce 



lim — 

_R^+oo R 



Ui{uj,-y'i,0)dy[. 



(4.12) 



This limit exists almost surely by the ergodic theorem, and the measurable function U[' 
satisfies Uj^ o Ty-^ = for all yi £ R. More precisely, one has 



lim — 

R-^oo R 



Ui{uj,yi - y[,0)dy[ 



29 



uniformly locally in yi. Indeed, for all |yi| < M, R large enough, 



R 



- (?7Kw,yi-y'i,0)-[/Ku;,-y'i,0)) dy[ 

I'M Y f-R+M 

<-/ \Ui{uj,-y[,Q\dy[ + - \Ui{uj,-y[,Q\dy[ 

-ft J-M -ft JR-M 



\Ui{uj,-y[,{)\''dy[ 



< C{M) 



1 



2R,2R) 
1 



1/2 



^i? \R>1 R Ji-R,R) J 

this last quantity vanishing uniformly as R +oo, almost surely. Then, to show that 
^«°2 = 0' integrate div Ui = over W^iuj, R) n {yi < 0}: 







which gives, as ii ^ +oo. 



div Ui 

■R~iLu,R)n{yi<0} 



Ui ■ n, 



d{TZ7{uj,R)n{yi<0}) 



Ur,= hm 



1 



R^+oo R 



Ui^i dy2 

r7{-R)n{y2<0} 



Ui^i dy2 

<T(0)n{y2<0} , 



If Uj^ were not zero, one would get from the left hand side that 
/ >C [ > C r dr ( [ |C/,,i| + / \Ui,,\] > SR^ 

JtZ{u),R) Jn-{LO,R) JRo \J(Ti-r)n{y2<0} J cr{r)n{y2<0} J 



for some 6 > and R> Rq large enough. This would be in contradiction with H2.8|) b. so that 
We now show that almost surely. 



1 

hm - / Ui{io,y[,0)dy[. 



Let Uj^ denote this limit, which exists almost surely for the same reason as Uj^. The con- 
vergence also holds in L'^{P) because Jq Ui{u>,y[,0) dy[ € L'^{P): 

21 



hm E 



lim E 

R^+oo 



R 



- I Ui{u;,yi,0)dyi-Uriuj) 



2- 

UiiLJ,yi,0)dyi-Urico 

R 



0, 



0. 



In particular, we have 



lim E 

R^+oo 



lim E 

R-^+oo 



— I Ui{uj,yi,0)dyi 







2- 

Uiiu;,yi,0)dyi 

'R 



E[|c/rMp] 
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since the law of Ui is invariant by translation in the yi-direction. Now, we write 

r2R 21 



E[|f/rHI 



^1 = hm E 



hm E 



I Ui{uj,yi,0)dyi 







+ E 



1 

2R 

R 



R 



Uiiu;,yi,0)dy-i 



+ E 



1 

2R 



2R 2- 

Ui{uj,yi,0) dyi 



R 







2R 



Ui{uj,yi,0)dyi / Ui{u,yi,Q) dyi 



R 



In this expression, the first two terms tend to the same limit |E[|C/j°°(cj)p]. Hence we deduce 
that the third term converges to ^E[|C/j°°(a;)p]. Then we have 



2 lim E 

R-*+oo 



1 

J Ui{uj,yi,0)dyi J ^Ui{uj,yi,0) dyi 



2E[\Ui 



oo|2i 



2 hm E 



i?2 



2R 



R 



R 



Ui{uj,yi,0)dyi / Ui{uj,yi,0) dyi 



0. 



To obtain the convergence of Ui to Uj^ as y2 — > +oo, we write by lemma ll. Ill 

Ui{L0,yi,y2)-Ur{uj)= [ G{y[,y2) {Ui{uj,yi - y[,0) - Ur{uj)) dy[ 



Giy[, y2) ^7 r - 0) - Ur{uj)dz dy[ 

oUi Jo 

diGiy[,y2)y[ l^J^'UiiLu,yi-z,0)dz-Ur{ioyj dy[ 



Let e, M > 0. There exists R = R{M) > such that 

V|yi| < M, y\y[\ > R, 
We deduce that 



\ rUi{uj,yi-z,0)dz-Ur{uj) 
Vi Jo 



1 fVi 



y'i\>R 



Vi Jo 



diG{y[,y2)y[ - / Ui{oj,y^ - z,^) dz - Ur{i^) \ dy[ 



< e. 



< Ce 



Moreover, 



/ diG{y'i,y2)y'i \\ r Ui{uj,yi - z,o)dz -uriuj)] dy[ 

J\y[\<R yVl Jo J 

\Ui{u;,z,0)\dz + R\Ur{u^)\] [ \diG{y[,y2)\ dy[ 



< G 



R+M 



R-M 

< C{R,M)\ 



\y[\<R 



0. 



^2 y2^+co 

This proves the almost sure convergence result. The quadratic convergence is simpler and 
left to the reader. 
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4.4 Some more estimates 

We will also need in the sequel the following estimates. 
Proposition 12 We have for all 1/2 > and all a G N^, 

1 



sup- 
R>1 ^ 



R 



\dyUi{uj,yi,y2)\ dyi < 00 a.s. 



and for \a\ > 1, 



sup^E/ \d^Ui{-,yi,y2)\^dyi 

R>1 R J-R 



0. 



(4.13) 



(4.14) 



Note that we could also prove an almost sure version of (|4.14p . as in the previous proof, but 
we do not need it. For q = 0, ()4.14|) holds if one replaces Ui by Ui — Uj^, as we have seen 
before. 
Proof. 

We have 

d''Ui(u;,yi,y2) = [ d'^Giyi - y[,y2)Ui{io,y[,0) dy[ 



hence we have E[|9"C/i(-, 0, 2/2)!] < 00. The law of (?"[/;(•, yi, 2/2) is independent of yi. Thus 
we can apply the ergodic theorem, which yields ()4.13|) . We now prove 1)4.14(1 . We have 

d''Ui{uj,yi,y2) = J d''G{y[,y2m{uj,yi-y[,0)-Ur{uj))dy[ 

y[did''G{y[,y2) (^^J\ui{ij,yi - z,0) - Ur{uj)) dz^ dy[. 

Let M > 0. We cut the previous integral in two pieces: \y[\ < M, \y[\ > M, and apply the 
Cauchy-Schwarz inequality in each piece. Recalling that G is homogeneous of degree —1, we 
obtain: 



9°t/K2/i,y2)P<^S^ / \d,d"G{y[,y2)\ {Ui{u;,yi - z,0) - Ur{u;)) dz 
^2 ■^\yi\<^^ -^0 



dy[ 



+ 



C 



\y[did''G{y[,y2) 



1 pi 



Vi Jo 



{Ui{u;,yi-z,0)-Ur{u;))dz 



dy'i- 



We now take the expectation, and obtain, due to the invariance of the probability measure 
P with respect to r^^ , 



E[|5"C/K-,yi,y2)r]<- 



G 



\did''G{y[,y2)\E 



y'l 



{Uii;-z,0)-Undz 



dy'i 



+ / \y[d,d^G{y[,y2)\^ \ T (^z(., -z, 0) - C/^ ) 
yP J\y[\>M 2/1 Jo 



dz 



dy'i- 
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Hence the quantity in the l.h.s. is independent of yi. Thus we have 



1 

sup-E/ \d''Ui{u;,yi,y2)\^dyi< 

R ^ J-R 



C 



2+2|q| 
2/2 



E 



M 



iUiiu;,-z,0)-Undz 



C 

c cs 



1 fy'i 



Vi Jo 



'iUi{u;,-z,0)-Undz 



2+2\a\ ' 2\a\ 

y2 2/2 



provided M is chosen great enough so that E|^ /(fi(C//(-,-z,0) - U^)dz\ < 5 for all \y[\ > 
M. The result follows. 



5 Justification of Navier's law 

This section is divided into two parts. First we derive an approximation u^pp of based on 
the boundary layer analysis, and prove theorem^ Then using this approximation, we justify 
Navier's wall law and prove theorem |S1 



5.1 Approximation of 

The approximation ti^pp reads 



u^pp{uj,x) = u {uj,x) + eUi ( w, — , — ) +sUu [ to,—, 



e e J \ ^ e / (5.1 



+ en^'^(a;, x), 



It takes into account the boundary layer terms, and includes corrections to recover the correct 
boundary and flux conditions : this is the role of n^'*". 

More precisely, since Vi^u does not converges to when 7/2 tends to ±00, we subtract f/^°^ 
from Ui,^ in H2.6() . Then we deal with the boundary condition. The term e{Ui — Uj^) is non 
zero on the upper boundary of Q^. We truncate this term replacing it by e(C// — U^)lx2<i- 
Hence we have the correct boundary condition, but a jump appears on Si. So we introduce 
a function evi{x) such that vi{xi, 1) = — Ui{^, ^) on Si, vi = on Sq, with div f ; = 
in J7, and we set vi = outside Q. Notice that e{Ui — U^) = o(e), thus evi = o(e) (roughly), 
so this term will not interfere in our estimate. 

Now on the lower boundary, we have e{Ui — U^) = —eU^. So we will again introduce a 
counterflow eq in $7. Unlike evi, eci is not o(e), thus we have to choose for q a solution of 
the following equations in : 

' - uAci + {u^ • V)q + (q • V)^^ + Vgi=0 

div q = 

' qIso = U^^ei, q|si = 
1 

q • ei dx2 = 
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(recall that = 0). The solution of this system is just a combination of a Couette flow and 
a Poiseuille flow: q = (1 — 4x2 + 3x2)C/;°Jei, Vgi = Qi>Uf° . We extend q outside by setting 
Q = above Ei and q = C/;°° below Sq. 

We proceed similarly with e{Uu—U^), introducing Vu and c^. We have Cu = (Sx^— 2x2)?7^ 
in fl, Cu = above Si and = below Sq. Finally, we add a small Poiseuille flow eOuq 
to correct the value of the flux of u'^app- We obtain: 

ulpp = u'^ + ^{Ui(^^) -U^y^,^i + e(Uu[^)-U^y^,>Q + e{vi + Vu + ci + Cu)+eev!^ (5.2) 

where the coefficient is defined by 

The value of ^ is independent of xi, because div u^app — 0- Easy recombinations show that 
the expansion (|5.2|) is of type (|5.H) . with 

and li^ = Q + c„ - - (according to (jHH) and (|^ ). 



5.1.1 Construction of and 



We explain in this paragraph the construction of the flow vi. The construction of Vu is 
analogous. We have to solve the following problem: find vi G Hl^J^Vi) such that 

' t;;(xi, 0) = on So 
vi{xi,l) = Ur -Ui{^, I) on Si 
div = in 0. 



Proposition 13 This problem possesses a (non unique) solution vi such that 



sup 

R>i 



R 



< oo a.s. 



and 



R> 



when e ^ 0. 



Proof. 

We shall find a solution in the following form: vi = V'^il). In terms of ip, the bound- 
ary conditions can be rewritten as ditf^ = 821!^ = on Sq, 91^^(3^1,1) = —Ui2{^A) and 

1 ' c c 

d2ip{xi, 1) = f/;,i(^, -) — on Si. Up to a constant, this is equivalent to ^ = d2ip = 
on So, ^P{x^,l) = - Ui^^i^,^ dx[, d2i'{xi,l) = f//,i(f ,i) - Ul^^ on Si. We search a 
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solution which is polynomial in X2'- let iIj{xi^ X2) = a{xi)x2 + b{xi)x2 + c{xi)x2 + li(xi). The 
boundary conditions on Sq imply c(xi) = d{xi) = 0, and on Si we get 

3a + 26 = f/,,i{^i,-) -C/g. 

This system has of course a unique solution {a{xi),b{xi)), and to obtain the required esti- 
mates, it is sufficient to show that 



1 

sup- 

R>1 ^ 



Xl 



Ul,2i^,'-)dx[ 



sup — E 

R>1 -K 



UiM^,hdx', 



2 



e ' e y 



r 7-00 1 1 ^ 



00112 



\uiAf,l)-ur,\ 



H'i{-R,R) 



< 00 a.s., 
> 0. 



The estimate on the quantity \\Ui{^, ^) — t^/°°||^3(_fl /j) follows from proposition El for we 



have 



l\\2 



(recall that Uj^ = 0). Hence we only need to deal with ■^||/q'^ ^i^'^^'e'^ ^^i\\'l'2(-r r)- 



order to do so, we write: 



Ui^2{^,\)dx\ = e Ui,2{yiA)dy 







xi/e 



Ui,i{^,y2)dy2-e 

a(-^)n{y2<l} 



Ui,i{0,y2)dy2 

a{0)n{y2<h} 



since Ui is divergence free on TZi n {0 < yi < ^ , 2/2 < Now the pointwise estimate is 
easy, for we have 



E 



Ui,i{0,y2)dy2 
<T{0)n{3/2<M 



< 00, 



hence we can apply the pointwise ergodic theorem. Next we turn to the estimate in expecta- 
tion. We have 



l/e 



1/s 



Ui,2{^,l)dx[=e {Ui,^{^,y2)-U[^,)dy2-e / {Ui,^{0,y2) - U^,) dy2 
+ £ 



Ui,i{^,y2)dy2- e 
<^{-^)n{y2<o} 



Ui^i{0,y2)dy2 
'7(0)n{y2<0} 



Then the invariance of P yields 
1 



E- 



2R 



e \ ^'\ui,i{^,y2)-Uri)dyr 



E 



L2{-R,R) 



C'\ui,i{Q,y2)-Uf^,)dy, 
Jo 
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Therefore it is sufficient to deal with the integral of the l.h.s. with R = 1. Let 6 > 0. We 
have 



E 



< e 



L2(-l,l) 



1/e 



E 



-1 
1 



< e 



Ui^ii^,y2)-Ul^,\^dxi 
fj'^ j'^\Ui,i{yi,y2)-U^^,\Uyi 

fl/e 

\Ui,i-Uri? + e 5 

J M 



dy2 
dy2 
dy2 



< sE , 

'[-l.l]x[0,M] 

<Ce + 6, 

where we have used successively the Cauchy-Schwarz inequality, the invariance of P and the 
estimate ()2.9() to make the expectation of J^-^ \ Ui^i — Uj^]"^ dyi less than 6 if y2 is larger than 
some constant M. Letting e decrease to 0, we deduce that the l.h.s. above tends to 0. It 
remains to deal with the integrals on cr{^) n {y2 < 0} and a{0) n {y2 < 0}. As above, the 
expectations of the quadratic mean on {—R, R) of these two integrals are equal, and equal to: 



E- 



2R 



<^{-^)n{s/2<0} 



'Xl 



y2) dy2 



L'2(~R,R) 



E 



E 



Ui^i{0,y2)dy2 

^(0)n{y2<0} 



Ui,i{yi,y2)dy2 

o-(s/i)n{y2<0} 



dyi 



\Ui,i{yi,y2)\ dy2dyi 

J (7iyi)n{y2<0} 



Putting together the four previous estimates, we can conclude that 

2 



sup — E 

R n. 



Ul,2C-^,l)dx[ 







0. 



L'^(-R,R) 



This ends the proof of proposition 1131 
5.1.2 Proof of theorem m 

This paragraph is devoted to the estimates 1)2.10(1 . It is enough to derive the first inequality, 
as the second one follows directly from Poincare inequality. In fact, we will show that 



sup — E 

R>1 ^ 



n^(-,R) 



where u^pp is given by (|5.1|) - ()5.2j) . Indeed, one verifies easily from the previous estimates that 



sup — E 

R>1 -K 



n^{-,R) 
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We write w = - ulpp = - v^pp, = - v^pp - ^^pp 
v = 0, 



^app — ■ It satisfies p.2|) with 

/(w, x) = 0, X in Q, 

f{uj,x) = ( -1217(1) {I + ee),o], X in n%oj)\ fl, 



and 



e 



|So,i + 1^ d2ivi + Vu + Ci + Cu)\eo^i + ^£04'- 



Finally, using that 

div (n° + + ^^^)) = div {{Ur + + u^) = 0, 
where is given in (|4.2j) . we write 



Proceeding as in section |21 one has 



ri+l 



r)+l 



dR I \Vw\'^ < C( I dR (|Gr+eVr) + e / dR 

+ / (|G|^ + ^Vl^) + / |VH^' 

which yields roughly, after integration with respect to uj: 

E r^\R [ [ \Vw\A < c(e [ + + eE [ |<^p 

Jr] \Jn^(R) J \ Jn^{v+^) ^So{»7+l)USi(r)+l) 

+ E [ iVu-l^ 



\4>\' 



Assuming that 



E / + eVi' < S{e) (t? + 1) , 6{e) — > 0, 



E 



|<^|2 < Ce^{ii + 1), 

/So,i(»?+l) 

we end up with a reverse Gronwall inequality: 



(5.3) 



/ e— >0 



Firj) < 7]{e)e\rj + l) + CF'{rj), F{rj)=E dR I [ |V 

and conclude as in section 13 It thus remains to establish bounds The second inequality 

is obvious using propositions 1121 and 1131 To control G, we first notice that 



\v'(^v'\' < C{ri + 1) sup Ib1li4(f,.(fc,fc+i)) 

n^iV+l) \k\<v+2 



< C{r^ + l)\\v'\\l, \\Vv 

uloc 
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by standard Gagliardo-Nirenberg inequality. Using estimates 1)3. 19() . H3.21() . we deduce 

lu^^v^P < Ce^(?7 + 1). 



/ 



Then, with notations of (|5.1() we write 



{u;,x)-e{ur{io)-U^{u;)-uHu;,x)) = e((c/-(^,^ 



Xl X2-I 

e e 



In particular, the estimate of G involves 

^{w,r?+l) 



E 



7 



£ e 



dx duj 



< C(r/ + 1)E e 



l/e+hu{0) 
MO) 



Uiiuj,0,y2)-UU^) dy2] du; = oil), e ^ 0, 



with the same manipulations as in the proof of proposition 1131 All other terms lead to similar 
inequalities and are left to the reader. This ends the proof. 

5.2 Effective wall law 

It remains to connect the approximation we have built to the Navier's wall law. Following 
|15j . such connection can be seen at a formal level, as 



u [uj, X) + eUl \^uj, — , — 



Xl X2-I 

e e 



+ £U^{ijJ, x), 



For instance, we obtain formally: 



+ a2C//,i(^,y,0 



which yields after averaging the approximate boundary condition, 

uflso = ea/<92Ui|so- 

Similarly, 

^ilsi = -ea„92wf|si- 

We now justify this formal computation and prove theorem [5J Let v'^ the solution of 



(|2.11|l - (|2.12j) . It is given explicitly by 



[x2- {x2 - eai)- 



1 {1 - 2eai){l - 2eau) 
3 ~ 



2(1 - e(Q; + Q„))y ' "y ■^^■"'^ 3 1 _ £^Q,^ _)_ Q,^ 

Then, we use the identity: 

u^uj, x) - v^uj, x) = u^iuj, x) - ulppiuj, x) + ulpp{uj, x) - v^uj, x) = {u^ - (w, x) 

+ e(^u^{uj,x) +eu\uj,x) - U^{uj) - U^{uj)^ -v'{uj,x) 
+ e [co, ^, ^) - f/r M) + e (Uu L, ^, ) - U^iu; 



38 



By theorem @J we have 



1 



sup — E 



R>i R \yJn{R) 
Then, simple calculations show that the function 

w = e {u°iu;,x) +eu\oj,x) - U^{uj) - U^{oj)) - v^{uj,x) 
which is the solution of 

w ■ Vw — vAw + Vq = 0, x G ri, 
div w = 0, X £ Cl, 



wi = 0, u;2('^,-)|eo,i = 0, 

o-(xi) 

dw{uj, •) 



satisfies 



w{uj, •) -eai{uj) 
w{uj, •) + eauiuj) 

sup — E 
it 



dx2 
dx2 



:e%H(4[/;~H + 2C/-iH), X2 = 0, 



R>1 -r*. \Jn{R) 

Finally, we proceed as usual to get 



0{e' 



. ^ E ( / UiL,^,^)- Uriu^) ' dx] dco = o(l) 
1 R \Jn(R) V e e / J 



sup — . 

_R>1 



and similar estimate for the upper boundary layer. This concludes the proof of theorem |21 
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